Energy distributions from three-body decaying many-body resonances 
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We compute energy distributions of three particles emerging from decaying many-body resonances. 
We reproduce the measured energy distributions from decays of two archetypal states chosen as the 
lowest 0+ and l+-resonances in ^^C populated in /3-decays. These states are dominated by sequential, 
through the *Be ground state, and direct decays, respectively. These decay mechanisms are reflected 
in the "dynamic" evolution from small, cluster or shell-model states, to large distances, where the 
coordinate or momentum space continuum wavefunctions are accurately computed. 

PACS numbers: 21.45.-l-v, 31.15.Ja, 25.70.Ef 



Introduction. Energy and momentum conservation 
guarantees that two particles, emerging from decay of 
a given quantum state, appear with definite kinetic en- 
ergies inversely proportional to their masses. In three- 
body decays the available energy can be continuously 
distributed among the particles. Prominent classical ex- 
amples are a-emission and beta-decay, respectively. Sur- 
prisingly enough the decay of a many-body quantum 
system into three particles has not been well described 
microscopically although discussed phenomenologically 
for various systems. The process depends on the ini- 
tial state and the dynamic evolution or equivalently the 
decay mechanism. 

This problem of three-body decay is common to sev- 
eral subfields of physics. The invention of Dalitz plots 
was an early attempt to classify the decay mechanisms 
by use of intermediate two-body doorway states [l| . The 
underlying dynamics in particle physics may be described 
as quark rearrangements. Similar decays occur in anni- 
hilation of a proton-antiproton pair from a Coulomb-like 
orbit into three mesons l2l|. In molecular physics an ex- 
ample is decay of excited states of the Ha-molecule into 
three hydrogen atoms [3]. In nuclear physics there ex- 
ists a large number of three-body decaying systems of 
disparate structures and decay mechanisms, e.g. various 
excited states of ^He, ^Li, ^^c, i7j^g_ y^^^^ j^j^j 

more 

high-quality experimental data become available in all 
subfields 0, [1,1^, Q and quantitatively accurate mod- 
els are needed to extract and understand the underlying 
physics. 

The purpose of the present letter is to compute the 
energy distributions for three-body decaying excited nu- 
clear many-body resonances. We shall assume that the 
resonances are populated in /3-decays and consequently 
only an outgoing flux is present. For reactions an in- 
going flux is required to provide the population of the 
decaying wavefunction. Such a generalization is eas- 
ily achieved by allowing initial conditions different from 
that of a resonance wavefunction. In all cases the ma- 
jor difficulty is to compute accurately the asymptotic 
large-distance three-body wavefunctions corresponding 



to genuine many-body resonances, which possibly dif- 
fer at small distances from cluster states formed by the 
emerging three particles [1] . 

At least four problems must be solved to overcome the 
difficulties, i.e. (i) the complex scaled resonance wave- 
functions must be accurately determined even though 
they vanish exponentially at large distances, (ii) the 
wavefunctions must be traced as they "evolve dynami- 
cally" from relatively small to asymptotically large dis- 
tances, (iii) the Coulomb problem of coupling continuum 
states at infinitely large distances must be solved, (iv) 
the mixture of two- and three-body asymptotics must be 
accurately determined. 

Theoretical framework. We use the hyperspherical 
adiabatic expansion method of the Faddeev equations 
combined with complex scaling |Tl| . The hyperradius 
p is the most important of the coordinates. For three 
identical particles of mass ttIq the definition is 
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where r,; is the coordinate of particle number i, R is the 
three-body center-of-mass coordinate and m^v = TOq/4. 

The two-body interactions are chosen to reproduce the 
available low-energy scattering data. A three-body po- 
tential with a range corresponding roughly to the three 
touching constituent particles is adjusted to reproduce 
the energy position of the many-body resonance under 
investigation. The wavefunction and the complex energy 
of the resonance are then defined. Coordinate and mo- 
mentum space angular wavefunctions are identical for a 
given total energy and an asymptotically large value of 
p [13 . Finally, the energy (Ea) distribution of the parti- 
cle is obtained by integrating the square of the resonance 
wavefunction over unobserved momenta. 

The numerical computations must first provide accu- 
rate wavefunctions from small to intermediate distances, 
where the relatively fast changes due to the crucial 
short-range interaction are completed. This is efficiently 
achieved with the Faddeev decomposition, and a hyper- 
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spherical harmonics basis size individually adjusted to 
the accuracy needed for the different partial waves. The 
smoother variation from intermediate to asymptotic dis- 
tances is analytical for short-range interactions The 
long-range Coulomb potentials can be treated numeri- 
cally precisely as any other potential in the transitions 
from small to intermediate distances. 
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FIG. 1: The real parts of the four lowest adiabatic effective 
potentials, including the three-body potentials, as functions 
of p for the 0"*" and 1"*" resonances of ^^C. The two-body in- 
teraction, obtained from FL9||. is a slightly modified version 
of the al-potential of [20]. The parameters of the three- 
body Gaussian potentials, S exp(— are 6 = 6 fm and 
-S = 20,92 MeV for 0+ and I+, respectively HH. 



However, special treatment is required in the variations 
from intermediate to large distances. The necessary ba- 
sis sizes become insurmountably large. Our solution is 
to compute an accurate wavefunction at intermediate, 
but relatively large, distances. This is achieved when an 
even larger hyperradius, compensated by a larger basis, 
leads to the same observables derived from the wavefunc- 
tions. This stability condition is difficult to reach when 
both three-body background continuum states are pop- 
ulated simultaneously with resonances in one or more of 
the two-body subsystems. At sufficiently large distances 
we can precisely identify these structures as components 
in the complex scaled wavefunctions related to different 
adiabatic potentials e.g. sequential decay proceeds 
through a potential approaching the corresponding com- 
plex two-body energy E2b whereas no intermediate 
structure is present for direct decay to the continuum. 

When the two-body intermediate states have large 
widths the related radial wavefunction decreases quickly, 
because then the adiabatic couplings are large. These 



states then dissipate fast into the continuum described as 
direct decay and the distinction becomes artificial. This 
process eventually happens for all sequential decays since 
the intermediate states are unstable. Classification into 
sequential and direct decay is related to the use of dif- 
ferent complete basis states, i.e. either two-body reso- 
nances and the third particle in the continuum or three- 
body continuum states. Interpretation as sequential or 
direct is then meaningful when a few states in one basis 
are sufficient while many are needed in the other. Such a 
distinction between paths producing the same observable 
is not possible in quantum mechanics. 

Illuminating archetypes. Accurate measurements of 
a-particle energy distributions are available from decays 
of 0+ and l+-resonances of [|, 0, [H, [11. The low- 
est 0+-resonance is often described as a cluster state 
[fl [Tsl [l6l . [TtI . [isj whereas the l+-resonance in contrast 
is referred to as a shell-model state without any signifi- 
cant cluster structure |JT, .15, ld[l3j[ll|- Furthermore, 
the decay mechanisms are known to be different @, [3]. 
These cases are therefore ideally suited as illustrations of 
the present novel technique. 
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FIG. 2: The a-particle energy distribution for the 1"*" reso- 
nance of ^^C at 5.43 MeV above threshold at an excitation of 
12.71 MeV. The energy is measured in units of the maximum 
possible, i.e. 2 x 5.43/3 MeV. The interactions are given in 
fig- [D The thick solid and dashed curves are for coordinate 
space wavefunctions at p = 70, 100 fm. The thin curves are 
contributions from separate adiabatic potentials. The mo- 
mentum space computation (described below) fall on top of 
the p = 100 fm curve (thick dashed). The histogram is the 
experimental distribution [l3l ]. 



In fig.[l]we show the lowest potentials where the attrac- 
tive pockets at small distance support the resonances and 
provide the small distance boundary conditions. As the 
hyperradius increases beyond the barriers, the potentials 
all decrease as 1/p due to the Coulomb repulsion. The 
structure at large distances is necessarily of three-body 
character since this is the boundary condition imposed by 
the measurement. In contrast, at small distances these 
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clusters overlap and the detailed description must use the 
nucleon degrees of freedom. The first adiabatic potential 
corresponds to the ^Be(0+) state and therefore associated 
with this sequential decay. We shall explore the, perhaps 
surprising, conjecture that the decay can be described 
almost entirely within the present cluster model. 

The 1'^ -resonance. With the potentials in fig. [T] wc 
show the computed energy distribution in fig. [2] for the 
1+ -resonance where sequential decay via the ^Be ground- 
state is forbidden. The asymptotic behavior is reached 
already for hyperradii larger than about 60 fm. The small 
variation of the distribution from 70 fm to 100 fm show 
the convergence and the stability. Higher accuracy is ob- 
tained at these distances with a moderate basis size than 
at larger distances where the basis quickly becomes insuf- 
ficient. Two interfering adiabatic potentials are necessary 
to reach the impressive agreement with the measured dis- 
tributions. It is remarkable that the cluster model pro- 
vides this accuracy in spite of the fact that the initial 
decaying state is a many-body resonance without any 
three-body structure. 




FIG. 3: The four lowest radial wavefunctions as functions 
of p for each of the two 0"'" -resonances of '^^C at 0.38 MeV 
and 3.95 MeV above threshold or at excitation energies of 
7.63 MeV and 11.2 MeV Ratios and the small-distance 
dominating wavefunctions are given for both states. 

To test the reliability we Fourier transformed the wave- 
function in two ways, i.e. first numerically with coordi- 
nates from p = to 100 fm and secondly by use of the an- 
alytic solution obtained from the parametrized adiabatic 
potentials which asymptotically are sums of l/p and l/p^ 
terms. The results are remarkably similar distributions 
and the analytic result is in fact indistinguishable from 
the curve for p = 100 fm in fig. [2l Small deviations from 
the measurements could be due to two-body interactions 
with resonance properties deviating slightly from the val- 
ues measured for ^Be. However, most of the differences 
are more likely due to uncertainties arising from accep- 
tance of the detectors used in the experiment. 



It is amusing to estimate that the sequential decay via 
the ^Be 2''"-resonance would produce a similar central 
peak of about twice the width. To reproduce the data 
strong interference would then be necessary. Good phe- 
nomenological reproduction of the data is obtained by 
R-matrix theory where the smaller width is explained 
due to preferentially populating the low energy tail of 
the *Be 2+ -resonance, and where effects of interference 
also play an important role 0. 
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FIG. 4: The a-particle energy distribution for the second 0"*" 
resonance of ^^C at 4.3 MeV above threshold at an excita- 
tion of 11.2 MeV 0, [l3|. The interactions in fig. [T]give an 
energy of 3.95 MeV above threshold We use 2.8 MeV 

to account for interference and beta-feeding distortion. The 
maximum energy for the most likely resonance position is then 
2 X 2.8/3 « 1.8 MeV. The histograms are experimental direct 
(the small part), sequential and total distribution. 

The 0^ -resonances. The complex scaled radial wave- 
functions are shown in fig. [3] for the two lowest O"*"- 
resonances. The largest probability is found at small 
distance and they all vanish with increasing hyperra- 
dius. Their relative sizes are fairly insensitive to vari- 
ations of the hyperradius at large distances where the 
energy distributions are determined. The first resonance 
is described by the first adiabatic component for all dis- 
tances whereas the second resonance changes character 
from small to large p from 4% to 75% of the first adiabatic 
component. This component, which dominates for both 
resonances at large distance, approaches the configura- 
tion of the O^-resonance in *Be with the third a-particle 
far away, i.e. s- waves between each pair of a-particles 
for each of the three Faddeev components. The result 
is an energy distribution with characteristic features of 
sequential three-a decay, i.e. a narrow high-energy peak 
and a distribution around one quarter of the maximum 
of the ^^C resonance energy. Unfortunately these com- 
puted distributions are not accurate because the two- 
body asymptotic behavior in these cases are not reached 
for p less than 100 fm. However, the method provides the 
amount of sequential decay and we can substitute the in- 
accurate component by the known two-body asymptotic 
behavior. The energy distribution from decay of the first 
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0+-resonance at 0.38 MeV is then seen from fig. [3] to be 
almost entirely determined by the first potential which 
means sequential decay. The direct decay is only about 
1% in agreement with the experimental upper limit p^ . 
This energy distribution is then in complete agreement 
with experimental data. 

The second 0+-resonance is also dominated by the first 
adiabatic potential at large distance. This is in striking 
contrast to the domination by the second potential at 
small distance. This is an example of the importance of 
the dynamical evolution from small to large distances. 
The result is about 75% sequential (first potential) and 
25% direct decay described by the other adiabatic poten- 
tials. In comparison with measurements complications 
arise for two reasons both related to the large width of 
the order 1 MeV. First, effects of energy-dependent feed- 
ing in the beta-decay populating the decaying state are 
substantial in the data fn, ■ Higher beta-energies are 
rather strongly favored resulting in distributions moving 
towards lower energies. Second, the experimental analy- 
sis is hampered by possible effects from other resonances. 
Their contributions are possibly not fully disentangled. 

The peak energy corresponding to the resonance po- 
sition is at about 2.8 MeV in the beta-feeding process 
0, [13] . We illustrate in fig. [4] the sequential part of the 
energy distribution of the first emitted a-particle by us- 
ing the Breit-Wigner distribution defined with the most 
probable position at 2.8 MeV and a width equal to the 
sum of the widths of the three-body decaying resonance 
and the intermediate two-body resonance. The two a- 
particles following from decay of *Be are uniquely related 
by kinematic conditions resulting in a peak at lower en- 
ergy. The large width of the three-body decaying reso- 
nance smears out the latter distribution. Between these 
two peaks appears the contribution of about 25% from 
direct decay described by the other adiabatic potentials. 
The inaccuracies in the computed distributions are first 
that deviations from the Breit-Wigner shape become im- 
portant for the large width of 1 MeV, and second that 
the fraction of sequential decay may be underestimated 
by perhaps 10% due to missing higher partial waves. 



In any case, the shape of the sequential decay via the 
^Be ground state is derived by precisely the same kine- 
matic conditions in both the computation and the ex- 
perimental analysis. The largest differences between the- 
ory and experiment is simulated by the shift of resonance 
peak energy. The agreement is rather good and only pos- 
sible due to the computed decay mechanism of dynamic 
evolution with hyperradius. 

Summary and conclusions. We have computed the 
energy distributions for three-body decaying many-body 
resonances. Combinations of short-range and repulsive 
Coulomb interactions are allowed. We conjecture, and 
show in specific cases, that the energy distributions of 
the decay fragments are insensitive to the short-distance 
many-body structure, but accessible in a three-body clus- 
ter model. The resonance structures may be completely 
different at small and large distances. This dynamic evo- 
lution is decisive for the decay mechanism. We separate 
components with two- and three-body asymptotics cor- 
responding to sequential and direct decays. This distinc- 
tion is crucial to obtain accurate wavefunctions at large 
distances. We test the method by comparing results from 
coordinate and momentum space. 

We illustrate by application to the archetypes of a- 
decaying 0"*" and l+-states in ^^C. The l+-resonance can- 
not be described as a three-body state but its decay 
proceeds directly into the three-body continuum. The 
two 0"'"-resonances both have substantial, but very dif- 
ferent, cluster components at small distances. However, 
they both decay preferentially through the same large 
distance structure best described as the O^-resonance of 
®Be. These sequential decays imply a total rearrange- 
ment of the second of these resonances from small to 
large distances. The accurately measured a-particle en- 
ergy distributions for all three resonances populated in /3- 
decay are reproduced remarkably well. Thus the method 
has past very severe tests. It is reliable and with predic- 
tive powers. 
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